Dynamics of spherical particles on a surface: 
About collision induced sliding and other effects 
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We present a model for the motion of hard spherical particles on a two dimensional surface. The 
model includes both the interaction between the particles via collisions, as well as the interaction 
of the particles with the substrate. We analyze in details the effects of sliding and rolling friction, 
that are usually overlooked. It is found that the properties of this particulate system are influenced 
significantly by the substrate-particle interactions. In particular, sliding of the particles relative to 
the substrate after a collision leads to considerable energy loss for common experimental conditions. 
The presented results provide a basis, that can be used to realistically model the dynamical prop- 
erties of the system, and provide further insight into density fluctuations and related phenomena of 
clustering, structure formations, and inelastic collapse. 

PACS numbers: 46.10.+Z, 83.I0.Pp, 46.30. Pa, 81.40.Pg 



I. INTRODUCTION 

In this paper we address the problem of the motion of 
a set of hard spherical particles on an inclined, in gen- 
eral dynamic surface. While there have been substantial 
efforts to understand in more details the problem of the 
nature of interaction of a single particle with the sub- 
strate [fl| — flOf| , these efforts have not being extended to 
the multiparticle situation. On the other hand, there has 
been recently a lot of interest in one jll] [h| or two 1 17 2^1 
dimensional granular systems. These systems are of con- 
siderable importance, since they provide useful insight 
into more complicated systems arising in industrial ap- 
plications, and also because of many fascinating effects 
that occur in simple experimental settings and theoreti- 
cal models. Theoretical and computational efforts have 
lead to results including density fluctuations, clustering, 
and inelastic collapse Jl^- p^ , ^6| -[2r| . Further, a system 
of hard particles energized by either an oscillating side 
wall, or by oscillating surface itself, has been explored 
recently p^ , ^6||20[| . This system, due to its similarity to 
one or two dimensional gas, appears to be a good can- 
didate for modeling using continuous hydrodynamic ap- 
proach 

It is of great interest to connect theoretical and com- 
putational results with experimental ones. Very recently, 
it has been observed experimentally that many complex 
phenomena occur in the seemingly simple system of hard 
particles rollin g and/ or sliding on a substrate. In particu- 
lar, clustering Qj^], and friction-based segregation 0] 
have been observed. While some of the experimental 
results (e.g. clustering) could be related to the theoreti- 
cal results |^0), there are still considerable discrepancies. 
Theoretically, it has been found that the coefficient of 
restitution, measuring the elasticity of particle-particle 
collisions, is the important parameter of the problem, 
governing the dynamical properties of the system. While 
the coefficient of restitution is definitely an important 
quantity, a realistic description of an experimental sys- 
tem cannot be based just on this simple parameter. As 



pointed out in [pl| , the rotational motion of the particles 
and the interaction with the substrate introduce an ad- 
ditional set of parameters (e.g. the coefficients of rolling 
and sliding friction), that have not been included in the 
theoretical descriptions of the system. 

Our goal is to bridge this gap between experiment 
and theory, and formulate a model that includes both 
particle-particle and particle-substrate interactions, al- 
lowing for a comparison between experimental and the- 
oretical results. Specifically, we address the phenomena 
of rolling friction and sliding, that lead to the loss of me- 
chanical energy and of linear and angular momentum of 
the particles. In order to provide a better understand- 
ing of the importance of various particle and substrate 
properties that define the system (e.g. rolling friction, 
sliding, the inertial properties of the particles), we con- 
centrate part of the discussion on monodisperse, hard 
(steel), perfectly spherical solid particles, moving on a 
hard (aluminum, copper) substrate pl|-p3|. However, 
through most of the presentation, the discussion is kept 
as general as possible, and could be applied to many 
other physical systems. Specifically, the extension of the 
model to more complicated systems and geometries is 
of importance, since most of the granular experiments 
involve some kind of interaction of the particles with 
(static or dynamic) walls. In particular, the discussion 
presented here is relevant to wall shearing experiments, 
where particle- wall interactions are of major importance 
in determining the properties of the system (see ]2q| and 
the references therein). 

In section || we explore all of the forces that act on 
the particle ensemble on a moving, inclined substrate. 
First we explore particle-particle interactions, and for- 
mulate a model incorporating the fact that the particles 
roll on a surface, and have their rotational degrees of 
freedom considerably modified, compared to "free" par- 
ticles. Further, we include the interaction of the parti- 
cles with the substrate, paying attention to the problem 
of rolling friction and sliding. The analysis is extended 
to the situation where the substrate itself is moving with 
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the prescribed velocity and acceleration. Because of the 
complexity of the interactions that the particles experi- 
ence, we first consider the problem of particles moving 
without sliding, and include the sliding at the end of the 
section. In section [II , we give the equations of motion for 



a particle that experiences collisions with other particles, 
as well as interaction with the substrate. In Section IV 
we apply these equations of motion to the simple case of 
particles moving in one direction only. It is found that 
many interesting effects could be observed in this simple 
geometry. In particular, we explore the effect of sliding 
both during and after the collisions, and give estimates 
for the experimental conditions that lead to sliding. Fi- 
nally, we give the results for the time the particles slide 
after a collision, as well as for the sliding distance, and 
for the loss of the translational kinetic energy and linear 
momentum of the particles. 



II. FORCES ON PARTICLES 

Particles moving on an inclined hard surface experience 
three kind of forces: 

• Body forces (gravity); 

• Forces due to collisions with other particles and 
walls; 

• Forces due to interaction with the substrate. 

In what follows we analyze each of these forces, with 
emphasis on understanding the interaction between sub- 
strate and the particles. While the analysis is kept as 
general as possible, some approximations appropriate to 
the problem in question are utilized, in order to keep 
the discussion tractable. In particular, the coefficient 
of rolling friction is assumed much smaller then (static 
and kinematic) coefficients of sliding friction. Further, 
in this section the particles are confined to move on the 
substrate without jumping; the experimental conditions 
under which this extra degree of freedom is introduced 
are discussed in section [V . Throughout most of this sec- 
tion it is assumed that the particle s are moving on the 
substrate without sliding; in section II C 3 we explain the 
conditions for sliding to occur. 

In order to formulate a model that can be used for ef- 
ficient molecular dynamics simulations, we choose rather 
simple models for the interactions between the particles 
and between the particles and the substrate. In modeling 
collisions between particles,, we neglect static friction, as 
it is often done f2q-|3lj|. On the other hand, the static 
friction between the particles and the substrate is of ma- 
jor importance, since it leads to rolling particle motion; 
consequently, it is included in the model. 



A. Body forces 

Here we consider only gravitational force that acts on 
the center of mass of the particles. It is assumed that 
there are no other (e.g. electrostatic) long range forces. 
In the coordinate frame that is used throughout (see Fig. 
1), the acceleration of a particle due to gravity, g, is 



»g = -j.gsin(6>) , 
where 9 is the inclination angle. 



(1) 




FIG. 1. The coordinate frame used in the paper. 



B. Collision forces 

There are many approaches to mode ling collision inter- 
actions between particles (see, e.g. [^,^2[ and references 
therein) . We note that rather complex models have been 
developed but choose to present a rather sim- 

ple one, which, while necessary incomplete, still models 
realistically collision between particles moving with mod- 
erate speeds. In the context of the particles moving on 
a substrate, it is important to realize that, even though 
particles are confined to move on a 2D surface, the 3D na- 
ture of the particles is of importance. Even if one assumes 
that the particles roll on a substrate without sliding, only 
two components (x and y) of their angular velocity, fi, 
are determined by this constraint. The particles could 
still rotate with O z , which could be produced by colli- 
sions (hereafter we use ft to denote the components of 
angular velocity in x — y plane only). We will see that 
rotations of the particles influence the nature of their in- 
teraction, as well as the interaction with the substrate. 

Normal force. Using a simple harmonic spring 
model |28| , p9| , ^2| , the normal force on particle i, due to 
the collision with particle j, is given by 



= [k(d ~ r id ) - 7j V m(v ij .,- • n)] n , 



(2) 



where k is a force constant, 
n = ri.j/nj, Vi.i = v, 



1,3 ~ I'ijli '!J ~ '! 3 ' 

ljj/.lj, . lj3 — y, - Vj, m is the reduced mass, 
and d = R{ + Rj, where Ri and Rj are the radii of the 
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particles i and j, respectively. In this paper, we assume 
monodisperse particles, so that d is equal to the diame- 
ter of a particle. The energy loss due to inelasticity of 
the collision is included by the damping constant, 7at. 
The damping is assumed to be proportional to the rel- 
ative velocity of the particles in the normal direction, 
n. While we use this simple linear model, the parame- 
ters k and are connected with material properties of 
the particles using nonlinear (Hertzian) model (see Ap- 
pendix |b]). We note that 7^ is connected with the coef- 
ficient of restitution, e„, by e„ = exp(— 7 n £ co ;/2). Here 
tcoi is the collision time and is approximately given by 
t C oi ~ ny/m/ (2k) (see Appendix |Aj). While more realis- 
tic nonlinear models lead to the velocity dependent co- 
efficient of restitution Q,^l|,^2), the linear model is sat- 
isfactory for our purposes, since we are interested in the 
collisions characterized by moderate impact velocities. In 
the case of particle- wall collisions, d — > d/2, fh — > m, and 



where r„ 



is the coordinate of the 



' *wall 

contact point at the wall. 

Tangential force in x — y plane. The motion of the 
particles in the tangential direction (perpendicular to the 
normal direction, in x — y plane), leads to a tangential 
(shear) force. This force opposes the motion of the inter- 
acting particles in the tangential direction, so that it acts 
in the direction that is opposite to the relative tangen- 
tial velocity, v* el) of the point of contact of the particles. 
Both translational motion of the center of mass, and rota- 
tions of the particles with component of angular velocity 



in the k direction contribute to 



thus 



u rel 



V 



+ 77 (n? + n?) 



(3) 



where § = (n • : 
particle i) by |2S 



-n • i). We model this force (on the 



r s = sign(-vl el )min (■y s m\v t rel \,v k \F L N \) 



(4) 



Here the Coulomb proportionality between normal and 
shear (tangential) stresses requires that the shear force, 
\F S \, is limited by the product of the coefficient of kinetic 
friction between the particles, z//c, and the normal force, 
F^|. The damping coefficient in the tangential (shear) 
direction, 75, is usually chosen as 75 = 7tv/2, so that 
the coefficients of restitution in the normal and shear di- 
rections are identical [p9| . An alternative method, where 
one models shear force by introducing a "spring" in the 
tangential direction, and calculates the force as being 
proportional to the extension of this spring, has been 



e.g. 



We neglect static fric- 



used as well (see, 

tion §§H|. 

The torque on the particle i, due to the force F|, is 
Ti = Xj xFj, where Xj is the vector from the center of the 
particle i to the point of contact, so Xj = —d/2 n. This 
torque produces the angular acceleration of the particle i 
in the k direction, flf = T*/I = —d/(2I) n x F|, where 
/ is the particles' moment of inertia. Recalling that F s 



is defined by Eq.(|]), one obtains (we drop subscript 
hereafter, if there is no possibility for confusion) 



n z = 



21 



(5) 



By direct integration, this result yields fl z . With this 
the relative velocity, -y* ei (Eq. (||)), is given, and hence 
we can calculate the tangential (shear) force, F s . 

Tangential force in the k direction. Since the particles 
are rolling, there is an additional force due to the rela- 
tive motion of the particles at the point of contact in the 
perpendicular, k, direction. Figure 2 gives a simple ex- 
ample of the collision of two particles with translational 
velocities in the i direction only. We model this force, 
FJj, which is due to rotations of the particles with angu- 
lar velocity, fi, in the same manner as the shear force, 
Fg. The force, F^, on the particle i due to a collision 
with the particle j, acts in the direction opposite to the k 
component of the relative velocity of the point of contact. 
Similarly to the "usual" shear force, we assume that the 
magnitude of F^j cannot be larger than the normal force 
times Coulomb coefficient; thus 



F C R = sign(-v* el )min(j s m\v* el \,i'k\F c N \)ic 



(6) 



and, for a general collision, v* el = + Flj) x k] • n. 

In the case of a central collision as shown in Fig. 2, v* el 
simplifies to 



u rel 



(7) 




FIG. 2. The collision between two particles, with the linear 
velocities in the i direction only. The direction of follows 
from the assumption that |v,-| > |vj|. The friction force f is 
explained in the following section. (The particles i and j are 
assumed to be in contact; for clarity reasons we show them 
separated). 



The torque acting on the particle i due to this force 
is T = x x Fjj. This torque produces an angular accel- 
eration ft = T/J. Assuming that there is no sliding of 
the particles with respect to the substrate, we obtain the 
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following result for the linear acceleration of the particle 



sign(-v* el )n 



(8) 



A 2D example given in Fig. 2 shows that of both 
particles, i and j, is in the direction of n, or — i direction. 

Let us also note that the force modifies the normal 
force, Fjv, with which the substrate acts on the particle. 
From the balance of forces in the k direction, it follows 
that the normal force is given by 



• N 



mg\iCOs(8) 



(9) 



The 



jump 



more details in section [V 



Fjj is discussed in 



condition, mgkcos(#) 

Here we assume nonzero F^v, 
and consider only the motion in x — y plane. 

To summarize, the collision interactions of the particle 
i with the particle j lead to the following expression for 
the total acceleration of particle i (in x — y plane) 



a ' — a N + a S + a R ) 



(10) 



where 








a N = 






m 




1 


a| = 


— i 




m 




R 2 


a R = 


I 



% = — min{^ s m\v z rel \,v k \F c N \)sign{-v z rel )n . 



Here b. c n is the acceleration due to the normal force given 
by Eq. (0), is the tangential acceleration due to the 
shear force, given by Eq. (Q), and a^ is the rotational 
acceleration due to the tangential force in the k direction, 
given by Eq. @. 



C. Interaction with the substrate 

The theory of rolling and sliding motion of a rigid body, 
even on a simple horizontal 2D substrate is complicated. 
For example, even though the question of rolling fric- 
tion was addressed long ago H , more recent works jj [Tc| l 
show that there are still many open questions about the 
origins of rolling friction; similar observation applies to 
sliding friction. In order to avoid confusion, we use term 
"friction" to refer to either static or kinematic (sliding) 
friction; rolling friction is considered separately. 

We approach this problem in several steps. After the 
introduction of the problem, we first consider a parti- 
cle rolling without sliding, with vanishing coefficient of 
rolling friction, fj, r . Next we present the generalization, 
that allows for nonzero /x r , as well as for the possibility 
of sliding. The substrate is assumed to move with its 
own prescribed velocity, vg, and acceleration, a$, which 
could be time dependent. The generalization to space 



dependent vs and as is straightforward, but it is not 
introduced for simplicity. Similarly, we assume that the 
substrate is horizontal; the generalization to an inclined 
substrate is obvious. 

Figure 3 shows the direction of the forces acting on 
a rolling particle. The friction force f , that causes the 
particle to roll, acts in such a direction to produce the 
torque, T, in the direction of the angular acceleration of 
the particle. Assuming that this friction force is applied 
to the instantaneous rotation axis, it does not lead to a 
loss of mechanical energy, as pointed out in Q. If fi r is 
zero, the particle will roll forever on a horizontal surface. 

On the other hand, the rolling friction force, f r , acts 
in such a way to oppose the rotations. Thus it produces 
the torque, T r , in the direction opposite to the angu- 
lar velocity of the particle. This torque could be under- 
stood if one assumes a small deformation of the substrate 
and/or particle, that modifies the direction of the rolling 
friction (reaction) force, f r , applied to the particle at a 
point slightly in front of the normal to the surface from 
the particle's center We note that this reaction 

force is actually our usual normal force, F^r. While we 
include the rolling friction in the discussion, we neglect 
the small modification of the normal force due to the 
effect of rolling friction. 



1. Rolling without sliding and without rolling friction 

In this work, we ignore the complex nature (see 
e.g. |lj|^l) of the friction force, and assume that there 
is a single contact point between a particle and the sub- 
strate, with the friction force, f , acting on the particle 
in the plane of the substrate, in the direction given by 
Newton's law. In order to calculate the acceleration of 
the particle, we use the simple method given in ]Ioj . The 
approach is outlined here, since in the later sections we 
will use the same idea in the more complicated settings. 

If the substrate itself is moving, the friction force 



fill 



is responsible for the momentum transfer from the sub- 
strate to the particle, where a is the particle acceleration. 
This force produces a torque (see Fig. 3) 



Rk xf = m 



(12) 



where fl is the angular acceleration of the particle. As- 
suming that there is no sliding, the velocity of the contact 
point is equal to the velocity of the substr ate, vg (this 
constraint will be relaxed in section [IC3, in order to 



model the more general case of rolling and/or sliding) 

v s = v + Rkxfl. (13) 

Multiplying Eq. ([l2|) by kx and using Newton law, we 
obtain 



4 



a = — - k x Ct . 

mR 



(14) 



Taking a time derivative of Eq. (p"3|), and combining with 
Eq. (Q), one obtains the following result for the acceler- 
ation of the center of the mass of the particle 



a = 



I 



(15) 



Since, for a solid spherical particle, I = 2/5 mi? 2 , we 
obtain a = 2/7 a^. So, the acceleration of a solid particle 
moving without rolling friction, or sliding, on a horizontal 
surface, is 2/7 of the acceleration of the surface, a,s [|10[. 




v,a 



FIG. 3. The forces and torques resulting from parti- 
cle-substrate interaction. The friction force, f, produces the 
torque, T, in the direction of the angular acceleration of the 
particle; the rolling friction force, f r , produces the torque, T r , 
in the opposite direction, so that it leads to the decrease of 
the particle angular velocity, f2. The deviation of f r from the 
k direction has been greatly exaggerated for the case of hard 
(e.g. metal) spherical particles. 



where the coefficient of rolling friction, // r (|v|), is de- 
fined by this equation, v = v/|v|, and F^r is the normal 
force. Alternatively, one could define the coefficient of 
rolling friction as the lever hand of the reaction force 
f r shown in Fig. 3 ||; for our purposes, the straight- 
forward definition, Eq. (|l6|), is more appropriate. For 
the case of steel spherical particles rolling on a cop- 
per substrate, the typical values of fi r are of the order 
of 1CP 3 |2l],|2^]. Realistic modeling of the experiments 
where rolling friction properties are of major importance 
(such as recent experiment p2| , which explores a system 
consisting of two kind of particles, distinguished by their 
rolling friction) requires accounting for velocity depen- 
dence of /i r = /i r (|v|). 

We note that there is an additional frictional force, 
which slows down the rotations of the particles around 
their vertical axes. While this additional force is to be 
included in general, we choose to neglect it here, since 
for the experimental situation in which we are inter- 
ested pl| , p2| ) the collisions between particles occur on 
a time scale that is much shorter than the time scale on 
which this rotational motion is considerably slowed down 
by the action of this frictional force (for other experi- 
mental systems, e.g. rubber spheres, this approximation 
would be unrealistic). 

To summarize, a particle rolling without sliding on a 
horizontal surface experiences two kind of forces: first, 
the surface transfers momentum to the particle, "pulling" 
it in the direction of its own motion and leading to the 
acceleration, a, given by Eq. (15). Second, due to the 
rolling friction, the particle is being slowed down, i.e. 
it is being accelerated with the acceleration, a#, in the 
direction opposite to the relative velocity of the particle 
and the substrate. 



3. Rolling with sliding 



2. Rolling without sliding with rolling friction 

The rolling friction leads to the additional force, re- 
sponsible for slowing down a particle on a surface. As 
already pointed out, this force produces the torque, T r , 
(see Fig. 3), in the direction opposite to the angular 
velocity of the particle, fi. The origins of this force 
are still being discussed. The effects such as surface de- 
fects, adhesion, electrostatic interaction, etc., that occur 
at the finite contact area between the particle and the 
substrate B, as well as viscous dissipation in the bulk of 
material have been shown to play a role. Fortu- 

nately, for our purposes, we do not have to understand 
the details of this force, except that it decreases the rel- 
ative velocity v = v - vj of the particle with respect 
to the substrate. The acceleration of the particle due to 
this force, a^, is given by 

a, = -^\)\F N \, (16) 
m 



Finally we are ready to address the problem of slid- 
ing. Sliding of a particle that is rolling on a substrate 
occurs when the magnitude of the friction force, result- 
ing from Eq. (|TT|), reaches its maximum allowed value 
|fmax|, where 

\fma X \ = M*|Fjv| ■ (17) 

Here Fjv is the normal force with which the substrate 
acts on a particle in the perpendicular, k, direction, and 
fi s is the coefficient of static friction between a particle 
and the substrate. Once the condition (|l7]) is satisfied, 
the friction force has to be modified, since now this force 
arises not from the static friction, but from the kinematic 
one. The direction of f is opposite to the relative (slip) 
velocity of the contact point of a particle and the sub- 
strate, u. Here u = u-vg, where u is the velocity of the 
contact point. The magnitude of f is equal to the prod- 
uct of the normal force and the coefficient of kinematic 
(sliding) friction, /Xfc 
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(18) 



The typical range of values of y, s and /ifc are 0.5 — 0.7 and 
0.1 — 0.2, respectively. In the subsequent analysis we ne- 
glect rolling friction, since the rolling friction coefficient, 
fj, r , is two orders of magnitude smaller than both /i s and 

The condition for sliding if there are no collisions. 
In this simple case, the friction force is given by 
Eqs. ( [TT[ [HI). From the the condition for sliding, |f| = 
\fmax\, where |f maa | is given by Eq. (]l7|), we obtain that 
the sliding occurs if 



as > 1 



mR 



I 

— * N\ 

m 



(19) 



As expected, if the substrate is accelerated with large 
acceleration, a particle slides. On a horizontal sur- 
face, the condition for sliding is |ag| > \{&s)min\, where 
|( a s)mm| = (1 + mR 2 /I) fi s g. For a solid steel sphere, 
fi s w 0.5, so |(as) mi „| w 1.75(7. We note that this result 
does not depend on the diameter of a particle. 



si = -y(k x f 



n x F%) , 



(21) 



where we concentrate only on the rotations in the x — y 
plane (the rotations characterized by il z enter into the 
definition of only). Following the same approach that 
led to Eq. @, one obtains the result for the linear ac- 
celeration 

_ [F c n + F g + (k - F|)n] + m(a s + a G ) 

TY18L — 775 . [ZZj 



1 



mR 2 



Similarly, one can solve for the friction force, f , that the 
substrate exerts on a particle 



N 



m(a G -a s )-2f(k.Fyn 



1 + 



I 



(23) 



For a solid particle, one obtains 



ma = — [F 



x t- F c s + (k ■ F^)n] + - m(a s + a G ) . (24) 



III. MOTION OF THE PARTICLES 

The preceding section gives the results for the forces 
that the particles experience, because of their collisions, 
as well as because of their interaction with the substrate. 
Now we consider the mutual interaction of these effects 
and give expressions that govern the motion of the par- 
ticles. 

Similarly to before, we consider first the case where the 
particles roll without sliding. Sliding is included in the 
second part of the section. 



Eqs. (|2^, effectively combine the acceleration due 
to the substrate motion, E q. (|l5| ), and the acceleration 
due to the collisions, Eq. (10). We note that, due to 



the interplay between angular and linear motion of the 
particles, the total acceleration, given by Eq. (p2|), is not 
simply the sum of the accelerations due to the collisions, 
gravity, and the frictional interaction with the substrate. 
These interactions are effectively coupled, and one should 
not consider them separately. 



B. Motion with sliding 



A. Motion without sliding 

In this section, we do not include rolling friction, since 
its effect is rather weak compared to the effects due to 
the collisions and the substrate motion. It is important 
to note that this approximation is valid only during the 
collisions; in between of the collisions, the rolling friction 
force has to be included, since it is the only active force 
other than gravity. 

The linear acceleration of a particle (in x — y plane) is 
given by 



'iv 



ma.Q 



(20) 



where F^ and F c s are the forces on a particle due to 
the collision, in the normal and tangential directions, re- 
spectively; f is the friction force, and a G is the acceler- 
ation due to gravity. Figure 2 shows a simple 2D exam- 
ple, where, for clarity, F c s , the rolling friction force, and 
the rotations, characterized by Q z , are not shown. The 
torque balance (generalization of Eq. (|l2|)) implies that 
the angular acceleration of a particle is given by 



The condition for sliding follows immediately from 



Eq. (§), 
yielding 



and from the sliding condition, |f| 



N 



I 



(k • Ffj)n + m(a G -a s ) 



N\ 



(25) 



m R- 
I 



If this condition is satisfied, then f is given by Eq. 

Let us concentrate for a moment on the condition for 
sliding of a solid particle, on a horizontal static substrate, 
and neglect F C R and F s . The sliding condition is now 
given by |F^r| = (1 + mR 2 /I) [i s \Fn\. As one would 
expect, this result resembles the condition for sliding due 
to the motion of the substrate, Eq. (p"9|), since the two 
considered situations are analogous (the acceleration of 
the surface as, plays the same role as the collision force, 
F^, scaled with the mass of a particle). Consequently, 
if the substrate is being accelerated in the direction of 
F^, a larger F^, is required to produce sliding. So, it is 
actually the relative acceleration of a particle with respect 



G 



to the substrate motion which is relevant in determining 
the condition for sliding. 

The effect of F^, on the sliding condition is more in- 
volved. Since Fjy is connected with F C R via Eq. (^), F C R 
modifies both sides of Eq. (p5|). The net effect of F C R is 
discussed in some details in section IV. 

If the sliding condition, Eq. d25)j 
to relax the no-slip condition, Eq. (1 
condition, we have 



is satisfied, one has 
Instead of no-slip 



u = v 



i?k x n 



(26) 



If u = vg, the sliding velocity, u = u — v s , vanishes. The 
equation for the "sliding acceleration" , u (relative to the 
acceleration of the substrate, as), follows similarly to 
Eq. (§), 



rau = 1 



f + F 



mR 2 



N 



(k • F^)n 



-m(a G - as) . 



(27) 



To summarize, the linear acceleration of a particle is 



given by Eq. (|20|), where the normal force, F^, is given 
by Eq. (§) and the tangential force, F|, by Eq. (f|). If 
there is no sliding, then the friction force, f, is given 
by Eq. ( |23| ); on the other hand, if the sliding condition, 
Eq. (p5|), is satisfied, f is given by Eq. Jig}). Further, a$ 
is the acceleration of the substrate, and a<2, the acceler- 
ation due to gravity, is given by Eq. (|l|). The angular 
acceleration of a particle, ft, is given by Eq. where 
the rotational force due to the collision, F R , is given by 
Eq. (ffl. Finally, the sliding acceleration, u, follows from 
Eq. (p7|). As mentioned earlier, the acceleration due to 
rolling friction, a^, is not important as long as much 
stronger collision, or friction forces, are present; it has to 
be added to the linear acceleration, a, for realistic model- 
ing of the motion of the particles between the collisions. 
The rotational motion of the particles, characterized by 
n z , enters into the model only by modifying the collision 
force between the particles in the tangential direction, 

F c s- 

The general expressions given in this section are used 
in the MD type simulations p3|, in order to simulate 
the motion of a set of particles on an inclined plane. In 
this paper, we apply the results to a simple setting, and 
obtain the analytic results which provide better insight 
into the relative importance of various interactions. This 
is the subject of the next section. 



IV. DISCUSSION 

The analysis of the preceding section gives rather gen- 
eral results, that provide all the information needed for 
modeling of the particles' motion. On the other hand, 
the complexity of the final results obscures simple physi- 
cal understanding. In this section, we concentrate on the 
particular case explored in recent experiments pUE2|, 



performed with steel particles on a metal substrate, and 
choose parameters appropriate to this situation. This 
system allows for significant simplifications, so that we 
are able to obtain rather simple analytic results. The as- 
sumptions which we use in what follows are summarized 
here for clarity: 

• Particles move just in one, i direction; 

• Particles are rolling without sliding prior to a col- 
lision; 

• The relative velocity of the particle prior to a col- 



lision, v. 







l v i ~ v ?l' i s assumed to be in the 
range 100 cm/s> v® el > 0.1 cm/s. It is assumed 
that the linear force model, Eq. (||), is appropriate 
for these velocities, but a nonlinear model is used 
to determine the approximate expressions for the 
force constant, k, and the damping parameter, 
(Appendices [a| and |b]) . For smaller relative veloci- 
ties, we will see that the interaction with the sub- 
strate substantially complicates the analysis. Still, 
this range of velocities is the most common one in 
the experiments |^l|,^2| , so we do not consider that 
this is a serious limitation; 

• The particles are assumed to be moving on a hori- 
zontal, static substrate; the analysis could be easily 
extended to other situations; 

• We neglect rolling friction, since its effect is negli- 
gible during a collision, or as long as the particles 
slide; 

• For simplicity of the presentation, we assume that 
the particles initially move with the velocities in 
the opposite directions; the final results are inde- 
pendent of this assumption. To avoid confusion, 
the particle i is always assumed to be initially ei- 
ther static, or moving in the — i direction. 

We are particularly interested in answering the follow- 
ing questions: 

• What is the condition for sliding to occur; 

• How long does a particle slide after a collision, and 
what is the distance traveled by a particle during 
this time; 

• How much of the translational energy and linear 
momentum of a particle are lost due to sliding. 

In order to fully understand the problem, we start with 
the simplest possible situation, and that is a symmetric 
collision of two particles moving with same speeds, but 
opposite velocities. Further, we assume the collision to 
be totally elastic. From this simple example, we conclude 
that the particle-substrate interaction is not important 
during a collision, at least not for the before mentioned 
range of particle velocities. Next we look into the case 
of more realistic, inelastic collision. Finally, we consider 
general inelastic, asymmetric collision of two particles. 
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A. Symmetric collisions 

Let us concentrate on the first part of a symmetric 
central collision of two particles, and analyze the forces 
acting on the particle i, as shown in Fig. 4a. The only 
collision force acting on the particle is F^, the substrate 
is static and horizontal, and the rolling friction can be 
neglected. In this case, Eqs. (Ell E3, \27\) simplify to 



ma = F 



ft 



mil = ( 1 



k x f 



mi?" 



f + F 



N 



Further, the friction force is given by 



(28) 
(29) 

(30) 



After understanding this basic situation, it is easier to 
understand the role of the remaining terms, given in sec- 
tion III B, but ignored in Eqs. ( p8f]3l| ). The contributions 
from gravity and the motion of the substrate are obvious. 
The analysis of the collision force in the tangential direc- 
tion, Fg, is similar to the one about F^, since the forces 
in the normal and tangential directions could be consid- 
ered independently. The contribution coming from F C R is 
discussed in the following sections. We note that in the 
case shown in Fig. 4 (the particles are initially moving 
with exactly opposite velocities), the contribution from 
F C R vanishes, since it is proportional to the relative ve- 
locity of the point of contact in the k direction. 



1. Symmetric elastic collision 



N 



f = 



1 



mR 2 
I 

-Hk\F N \\i 



if |f| < fi s \F N \ 
otherwise 



'N 



Let us assume that the sliding condition, f = /i s 
is satisfied. By inspecting Eq. ([5(j|), we observe that the 
first term on the right hand side is the one that decreases 
the sliding acceleration continuously, and possibly brings 
the particle back to pure rolling, since it acts always in 
the direction opposite to the sliding velocity. Because 
of the constraint on the friction force given by Eq. (31), 
and /ife < /is, the right hand side of Eq. (|30| ) gives a net 
contribution in +i direction. So, when sliding begins, 
the particle i experiences the sliding acceleration in the 
i direction, leading to the sliding velocity, u, in the same 
direction, as shown in Fig. 4a. In other words, for this 
situation, the particle i is still moving to the left, with 
angular velocity in the — j direction, but it is sliding to 
the right, with sliding velocity u. Let us also note that 
f slows down the angular motion of the particle, as it 
can be seen by inspection of Eq. (p9]). Next, we consider 
the typical situation during the second part of the colli- 
sion, when the particles are moving away from each other 
(Fig. 4b). Analysis shows that almost all of the conclu- 
sions about the situation depicted in Fig. 4a extend to 
this situation; in particular the directions of the sliding 
velocity, u, and the friction force, f , are the same. 




FIG. 4. The forces considered in this section (x — z plane). 
Part a) shows the first part of the collision, when the particles 
are still moving towards each other, and b) shows the second 
part of the collision. For clarity, only forces on the particle i 
are shown. The rotational motion is discussed in the text. 



Let us define the compression by i = (d - 



' 1,3 



)/d, 



(31) where r^j is the distance between the centers of colliding 
particles, and d is the particle diameter (x > is required 
if a collision occurs) . The maximum compression is given 
by (Appendix |a]) 



J rel 



duo ' 



(32) 



where, for a symmetric collision, the relative velocity 
v rei = and v° is the initial speed of the particles. 
For simplicity, we use scalar notation and the sign con- 
vention that +sign of the translational velocities refers to 
the motion in +i direction, and -(-sign of the angular ve- 
locities/accelerations to the rotations in +j direction. In 
obtaining Eq. (|32|), it was assumed that only the collision 
forces are important in determining x^ax- More careful 
analysis given in Appendix [S] provides justification for 
this assumption. The natural frequency, loq, associated 
with the linear force model specified by Eq. (|^), is given 
by Wq = 2k /to. It is related to the duration of the col- 
lision via t co i = tt/luq. The nonlinear force model (see 
Appendix 0), predicts that ojq very weakly depends on 



J rel ' 



we observe that typically ujq 



1CP J s, so that for 
x ss 2 x 1CP 4 . Figure 5a shows the 
on the relative initial velocity of the 



10 cm/s, : 
dependence of x^ 
particles. 

It is of interest to estimate the range of the compres 
sion depths for which the sliding condition, |f | = /i^F^I 
is satisfied, leading to sliding of the particles with re 
spect to the substrate 
(1 + mR 2 /I) [i s mg (see Eqs. (| 
sion for the normal force, Eq 



condition is satisfied for x^ 
Appendix ^|) 



Sliding occurs when |F^| > 
|3l])). Using the expres- 
), we obtain that this 




> x 



slip 



> T U 
_ x~ 



mR 2 \ /i s <7 



dJ 1 



where (see 



(33) 



Using the values of the parameters as specified in Ap- 
pendices [§] and |^, we note that for v® el sa 10 cm/s, 
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x min ~ 3 x 10 . For d — 4 mm, x^ min is on atomic 
length scale, so we conclude that a particle slides during 
almost all of the course of a symmetric, elastic collision. 
The results for x® min are shown in Fig. 5b. The depen- 
dence of u>o on v° (see Appendix f|) , leads to the increased 



values of x" min as v 
The compression, 



is reached at time t? 



measured from the beginning of the collision (see Ap- 
pendix ^|) 



r • =1 



mR \ ji s g 



(34) 



For the choice of parameters as given in Appendix H, we 
obtain t^ in » 10~ 8 , and t® min jt co i w 5 x 10~ 4 , confirm- 
ing our conclusion that sliding with the respect to the 
substrate is the dominant motion of the particles during 
a symmetric, elastic collision. 

From Eq. (^) we can also deduce under what condi- 
tions sliding occurs. Obviously, we require that x® min < 



Using Eqs. 



we obtain that the initial ve- 



locities of the particles have to satisfy v® e[ 3> v b , where 
(see Appendix y) 



iR 2 



Vs9 
w 



(35) 



For the set of parameters given in Appendix [b|, this ex- 
pression yields very small value, v b ss 10~ 2 cm/s. So, 
sliding occurs during almost all collision occurring in typ- 
ical experiments [|2l| , [22[ . 

Let us now look into the rotational motion of the parti- 
cles during a collision. The friction force is the only one 
which produces angular acceleration. Without loss of 
generality, we consider the particle i, which is assumed 
to move initially in — i direction, with angular velocity 
fi? = — v° /R. Integrating over the duration of the colli- 
sion gives the result for the angular velocity of the par- 
ticle at the end of the collision (at t = t co i) (see Ap- 
pendix |j) 



n 



fo_ 



mR 



— x 

mR 2 



Vs9 



1 



(36) 



The sliding velocity of the particle i at t = t co i, u{° = 



Rnj v (see Eq. (EQ)), now follows 



raR 2 g 



J rel 



TT^fc 



I U!() 

mR 2 



x 



1 



2u° 



t UJ 



1 



-/J, s — (Ik 



(37) 



This is important result, since the sliding velocity of a 
particle at the end of a collision determines the energy 
and momentum loss due to the slidingafter the collision. 
Using the parameters as in Appendix H, we note that the 



contribution of the second term in Eq. fl37p is approxi- 
mately 0.01 cm/s. So, we conclude that the frictional 
interaction of a particle with the substrate during a colli- 
sion only very weakly influences the sliding velocity of a 
particle at the end of the collision. Similarly, the angular 
velocity is only slightly modified, as shown in Fig. 4. A 
particle exits an elastic, symmetric collision with the an- 
gular velocity which is almost equal to its initial angular 
velocity, resulting in the sliding velocity equal to twice of 
its initial translational velocity. 

The fact that the particle-substrate interaction is neg- 
ligible during a collision follows also from a simple energy 
argument. Figure 5 shows that the maximum compres- 
sion depth is of the order of 10 -4 , in units of the par- 
ticles diameter. So, the order of magnitude of the ratio 
of the energies involved in the particle-substrate interac- 
tion, Ep- S « fi s mgdx max , and of the energy involved in 
the collision itself, E co u « k(dx max ) 2 /2, is given by 



E„ 



E n 



{J. s mgv° el tcoi _ gtcoi 
m(v° rel ) 2 



(38) 



u rel 



where Eq. ( |32] ) has been used. For v® el — 10 cm/s, we 
obtain E p _ s /E co i « 2 x 10~ 3 . Clearly, for all collisions 
characterized by very short collision times (equivalently, 
small maximum compression depths), this ratio is very 
small number, assuming common particle velocities. Cor- 
respondingly, the particle-substrate interaction during a 
collision influences very weakly the dynamics. Consider- 
able modification of this estimate could be expected in 
the case of "softer" collisions, where both the duration 
of a collision and maximum compression depth are much 
larger. 



2. Symmetric inelastic collision 

Inelasticity of a collision introduces damping param- 
eter, 7jv, which is related to the material constants in 
Appendices |X| and [b]. The damping is directly con- 
nected with the coefficient of restitution e n by jn — 
— 2/i co ; ln(e„) (Appendix [b]). The collisions of steel 
spheres are rather elastic (typically e„ « 0.9), so we 
are able to introduce a small parameter, e = ~/n/u>q s» 
— 2/tt(1 — e n ) <C 1. In what follows we perform consistent 
perturbation expansions of the equations of motion, and 
include only the corrections of the order O(e). For the 
completeness, we also include the terms due to the inter- 
action with the substrate, even though we have already 
shown that this interaction is not of importance for the 
physical situation we are interested in. 

The maximum compression is now given by (see Ap- 
pendix ^) 



dujQ 



(i-§« + °C^>) 



(39) 



Figure 5a shows the result for x max , for a few values of 
e n , using the parameters given in Appendix 0. In Ap- 







pendix |B] it is shown that the uiq and t co i are weakly de- 
pendent on the initial velocity (t co i ~ (i; ) -1 / 5 ), so that 
the results for the maximum compression scale with the 
initial velocity as x max ~ ( w rei) 4 ^ 5 ' resulting in the slight 
curvature of the x max curves in Fig. 5a. In Appendix ^ 
it is shown that for inelastic collisions, the time, t m i n , 
at which sliding starts goes to zero, since the corrections 
due to damping are typically stronger than corrections 
due to the particle-substrate interaction. Consequently, 
the angular acceleration, given by Eq. (p9|), is constant 
during whole course of the collision. For the particle 
i, Qi — fikgmR/I, so that at the end of the collision 

(t = tool) 



a J = o 



,, t mR fj, k gn _ ^ 
I w 



(40) 



where = —v°/R. The sliding velocity of the particle 
i at t = t co i is given by 



mR 2 ^fc<?7r 1 

~1 LJ~ Q 2 



(1 + e n )v° rel . (41) 



Similarly to the discussion following Eq. ( |37| ) , we observe 
that the friction during a collision leads to negligible cor- 
rections. In what follows, we ignore these corrections, 



and assume il 



S 



and it. 



(1 



«/2. 




e =1.0 



v (cm/s) 



v (cm/s) 



FIG. 5. a) The maximum compression, x max ; b) the min- 
imum compression, x^im required to produce sliding. Both 
quantities are scaled with a particle diameter. For inelastic 
collisions, x min = (see text). 



3. Sliding after a symmetric collision 

The particle i exits a symmetric collision with the 
translational velocity v{ — e n v° (in +i direction), and 
with the sliding velocity u{ , given by Eq. (fiTl). After 
the collision, it experiences the friction force, resulting 
in the sliding acceleration iii = —(1 + mR 2 / I)[ikg, as it 
follows from Eq. (|30|), where F^, is now absent. This 
friction force is present as long as the sliding velocity is 
nonzero. It slows down the particle, and leads to the cor- 
responding loss of the translational kinetic energy and 
linear momentum. Neglecting rolling friction, we obtain 
the result for the time, t s , measured from the end of the 



collision, when sliding stops (due to the symmetry this 
result is the same for both particles) 



\Ui\ 



f(l + e> r °e, 



1 



mR? ' 



(42) 



The translational velocity of the particle i at the time t s 
is given by 



v l {t = t s ) = 



1 



J rel 



1 I 



1 



1 



T.R 2 

I 



(43) 



The angular velocity of the particle i at this time is 
fif = Qi(t = t s ) = uf/R, since the particle does not slide 
anymore. We observe that the translational motion of 
the particles is considerably slowed down due to the fric- 
tion force; for solid spheres, and e„ = 0.9, \v%\ ~ 0.36 v°, 
(k = Equation (^) gives that for I/mR 2 > e„, 

vf is negative, meaning that the particle is moving back- 
wards at the time when sliding ceases. For almost elastic 
collisions of solid particles this condition is not satisfied, 
so the particles are still moving away from the collision 
point at time t s . 

Until the time t s , each of the particles travels the dis- 
tance s away from the point where the collision has taken 
place, given by 



8/i fc3 (1 + ^f) 



r(l+e„) 



1 



mR 2 
2^— e, 



(44) 



Figure 6 shows the results for t s and s. For v® el = 10 
cm/s, and e n = 0.9, the particles slide during the time 
t s w 0.03 s, and s « 0.1 cm. These results compare well 
with preliminary experiments. More precise analysis and 
the comparison with experimental results will be given 
elsewhere [E3F 
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FIG. 6. The time, t", until which a particle slides after the 
collision, and the distance, s, traveled during this time. Here 
v — Vg el is the initial relative velocity of the particles. The 
parameters are as specified in the text. 
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4- The loss of the translational energy and momentum due 
to sliding in a symmetric collision 

Let us define the energy loss due to sliding, AE s u p , as 
the difference between the translational kinetic energy of 
a particle just after it has undergone a collision, and its 
translational kinetic energy at time t s after the collision 
(scaled with the reduced mass) . Therefore, 



AE 



slip 



(45) 



The relative loss of energy is defined as AE s u p = 
AE s u p /E , where E = {v ) 2 . We note that very little 
energy is lost due to sliding while the collision is taking 
place. The sliding loss of energy occurs after the collision, 
and it is equal to the work done by the friction force. Us- 
ing Eq. (p3|), we obtain the result for the relative energy 
loss due to sliding, 



AE sHp = (1 + e„)- 



1 + 2e„ 



niR- 
I 



(1 



(46) 



Figure 7a shows AE s u p for a range of values of e„ (as- 
suming solid spheres). In the limit of an elastic collision, 
e„ — ► 1, and we obtain AE s Up ~ 0.8. So, a solid par- 
ticle loses approximately 80 % of its initial translational 
kinetic energy due to sliding, in a completely elastic sym- 
metric collision. 

It is also of interest to compare the relative energy loss 
due to sliding, AE s n p , with the energy loss due to inelas- 
ticity of a collision, AE co i. The latter is simply given by 
(1 — e^) (we neglect the small loss of energy due 



col 



AE, 

to interaction with substrate during a collision) , thus 



AE, 



col 



AE. 



= 1 



slip 



mR 2 



1 - e„ 



- 1 + 2e, 



mR 2 



(47) 



The result for AE co i/ AE s u p is shown in Fig. 7b. We 
observe that in the limit of low damping, the sliding is 
the main source of energy loss. This conclusion is in- 
dependent of the initial particle velocity or the particle 
diameter. 



b) 



0.4 ^ 




0.8 0.85 0.9 0.95 

e„ 




1.8 0.85 0.9 0.95 

e„ 



FIG. 7. a) The loss of the energy and momentum due 
to sliding, b) The ratio of the loss of the mechanical en- 
ergy and linear momentum due to inelasticity of a collision 
(AE co i, Apcoi), and due to sliding (AE s ii P , Ap 3 u p ). 



Similarly, the linear momentum lost due to sliding in 
a symmetric collision (relative to the initial momentum) 
is given by 



Ap s iip 



1 



mR? 
I 



(48) 



so (in an elastic symmetric collision), a solid particle loses 
approximately 60 % of its linear momentum because of 
sliding (Figure 7a). The ratio of the loss of the linear 
momentum due to inelasticity of the collision, defined by 
Ap co i = (1 - e„), and Ap slip , is given by 



Ap, 



col 



Ap s up 



= 1 



mR 2 \ 1 



1 + e r , 



(49) 



[We note that Ap co i is the loss of linear momentum of 
one particle in the lab frame; inelasticity of the collisions 
conserves the linear momentum of a pair of colliding par- 
ticles in the center of mass frame.] This result is shown 
in Fig. 7b. Similarly to the energy considerations, we ob- 
serve that for small damping, sliding is the main source 
of momentum loss. 



B. Asymmetric collisions 

Next we consider a central collision between par- 
ticles moving with different speeds, such as the one 
shown in Fig. 2. On a horizontal static substrate, the 
Eqs. (pii], |22], p7| ) now simplify to (index i emphasizes that 
the particle i is being considered) 



a = -|(kxf t + nxF^) , 
mm = (l+ — I fi + F c N/j 



(50) 
(51) 

^(k-F^)n. (52) 



Further, the friction force is given by 

-inf (k-F^)n 



Fc mR 
NA 



ft = 



-Mfc|Fjv,i|u, 



mR? 
I 



if |$ | < Ms|Fjv,i 
otherwise 



(53) 



The analy sis of a symmetric collision, given in sec- 



tion IV A , shows that the frictional interaction of the 
particles with the substrate during a collision can be ne- 
glected. We use this result in the following discussion 
and neglect ft in the analysis of the collision dynamics 
of an asymmetric collision. This frictional interaction is, 
of course, included in the analysis of the particles' mo- 
tion after a collision, since it is the only force acting on 
a particle on a static, horizontal substrate. 

Normal force is being modified due to F C R i , so that 



(mg-k-F^)k 



(54) 



11 



In Appendix |E]it is shown that, for typical experimental 
velocities, the corrections of F C R i due to its cutoff value 
(see Eq. (0)), could be ignored, since the cutoff leads 
to 0(e 2 ) corrections of the final angular velocity of the 
particle. So, we take i to be given by (see Eqs. (^|, [?])) 



7s 



R,i 



mR&tli + flj) j]k 



(55) 



during the whole course of a collision. The damping pa- 
rameter, 7s, is kept as a free parameter for generality 
(usually it is given a value 75 = 7at/2 p9|]). Only con- 
straint on 75 is that 75 /u>q <C 1, so that the coefficient 
of restitution is close to 1. 

The force F„ ,- modifies the rotational motion of the 
particle i. In Appendix M it is shown that the angular 
velocity of the particle at the end of a collision (t = t co i = 
tt/loq) is given by 



n{ = n? - c(n? + a!}) 



(56) 



where C = nmR 2 j s I '{2Jw ) = 0(e) < 1. Equation (|5|) 
is correct to the first order in e. Using this result, and 
the translational velocity of the particle i at t = t co i 
(Eq. (A7)), we obtain the sliding velocity of the particle 
i at the end of the collision 



1 



< =-o(l+^)(v?-v°) + C(v°+v«) 



(57) 



This result generalizes Eq. (ffl|), that gives the sliding 
velocity of the particles undergoing a symmetric colli- 
sion (the particle-substrate interaction during the colli- 
sion has been neglected). The tangential force, F^, leads 
to the last term in Eq. ( p7| ) , modifying the sliding velocity 
in an asymmetric collision. This modification depends on 
|v? + v!-|, which measures the degree of asymmetry in a 
collision. 

In order to exemplify the physical meaning of these re- 
sults, let us consider for a moment completely asymmet- 
ric case: a particle moving with initial velocity v° and un- 
dergoing elastic collision (7^ = 75 = 0) with the station- 
ary particle i. In this case, we obtain = 0, = — v°. 
So, the particle j is stationary immediately after the col- 
lision, but its rotation rate is unchanged (since in the 
limit 7s = 0, FJj vanishes, and the interaction with the 
substrate has been neglected), so that it has the sliding 
velocity equal to the negative of its initial velocity. Let 
us now consider the particle i. Its translational veloc- 
ity and sliding velocities are the same, ~v{ — u{ = v°, 
since immediately after the collision this particle has the 
translational velocity equal to the initial velocity of the 
particle j, but zero rotation rate. 

"Jumping" of the colliding particles. Let us finally 
address the assumption that the particles are bound to 
move on the substrate. From Eq. (|54|) we observe that, 
for large positive k • F R , this assumption could be vio- 
lated. The estimate is given in Appendix |], where it is 
indeed shown that a particle colliding with a slower par- 
ticle typically detaches from the substrate. Fortunately, 



the motion of a detached particle in the k direction is 
limited by very small jump heights, so that the modifi- 
cations of the results for the dynamics of the particles in 
x — y plane are negligible. On the other hand, the fact 
that a particle is not in the physical contact with the 
substrate during a collision simplifies the analysis of the 
collision dynamics, since particle-substrate interaction is 
not present. We note that we are not aware that detach- 
ment has been observed in the experiments performed 
with steel spheres moving with moderate speeds pl|p2| . 
Since this effect provides direct insight into a collision 
model, it would be of considerable interest to explore 
these predictions experimentally. 



1. Sliding after an asymmetric collision 

After a collision, the particles experience the friction 
force, which produces the sliding acceleration and mod- 
ifies the translational velocity. Figures 8-11 show the 
results for the time that the particles spend sliding, for 
the sliding distance, and for the changes in their transla- 
tional kinetic energy and linear momentum. All of these 
results depend only on the sum and difference of the ini- 
tial velocities of the particles. We define 



v m = (v? - v°) • i , 



(58) 



and show the dependence of our results on these two 
quantities. Since some of the approximations involving 
the rotational motion of the particles during collisions 
(see Appendix ^) are not valid in the limit \v m \ <C \Vp\, 
we do not consider the case \v m \ ~ (which occurs when 
the initial velocities of the particles are almost the same) . 
This is the only imposed restriction. 

Using Eqs. ([52], |57|), we obtain the time when sliding of 
the particle i stops(measured from the end of a collision) 
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(59) 



Figure 8 shows the result for the sliding time for fixed 
e n and C, as a function of v m and v p . For v p — 0, we 
retrieve the results for the symmetric collision, shown in 
Fig. 6. We observe that t\ just very weakly depends 
on v p ; this dependence disappears in the limit of zero 
tangential damping (C = 0), as can be seen directly from 
Eq. H). 

The translational particle velocity at t = if is = 



+ a»if, where a, = -^fegu 
), we obtain 
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During the time tf, the particle i translates for the dis- 
tance | Si | from the collision point, where Sj = (v{ + 
vf)£?/2. Figure 9 shows |sj|; contrary to the sliding time 
tf, the sliding distance does depend on the asymmetry 
of a collision. This dependence is present since |sj| is 
a function of both translational and sliding velocities of 
the particle i. On the other hand, tf depends only on the 
sliding velocity of the particle. 




FIG. 8. The sliding time t" of the particle i. The solid line 
shows the result for a symmetric collision. Here e„ = 0.9 and 
= 0.13. 

An interesting effect can be observed in Fig. 9: there is 
a particular combination of the initial particle velocities 
that gives vanishing sliding distance. The meaning of 
this result is that the particle returns to its initial position 
exactly at the time tf after the collision; this occurs when 
vf = -vf . Using Eqs. @ |A7|), we obtain the condition 
for zero sliding distance in terms of the initial velocities 
of the particles 
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(e„ + 2) +e„2f--2C 
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(61) 



For a completely elastic collision of solid particles, we ob- 
tain v? = — 6v°. Equation ( |6l| ) gives clear experimental 
prediction which can be used to explore how realistic the 
collision model is. 



2. The change of the translational kinetic energy and 
momentum due to sliding 

In this section we give the final results for the change 
of the translational energy and the linear momentum of 
the particles due to sliding after a collision. This results 
assume that the particles slide the whole distance s, so 
that there are no other collisions taking place while the 
particles travel this distance. Consequently, for a system 
consisting of many particles (as in pl[p2[ ), the change 
of the translational energy due to sliding depends on the 
distance traveled by the particles in between of the col- 
lisions, I. When I is on average much larger than the 



sliding distance, s, one could consider modeling the effect 
of sliding using "effective" coefficient of restitution , 
which we derive below. In this case, we find that this 
"effective" coefficient of restitution depends only on the 
usual restitution coefficient, e n , and on the geometric 
properties of the particles. On the other hand, if I » s, 
this "effective" coefficient of restitution will depend also 
on the local density and velocity of the particles. Wc 
explore this effect in more details in p4[ . 




FIG. 9. The sliding distance s of the particle i (|si| in the 
text). The solid line shows the result for a symmetric collision 
(e„ = 0.9 and C = 0.13). 

The change of the translational energy of the particle 



I, AEl Up , is defined as AEl Hp = | vf |^ - |v?|^. The trans- 
lational velocity of the particle when it stops sliding, vf, 
is given by Eq. (p0|), and the velocity of the particle at 
the end of collision, vf , is given by Eq. (|A7|). 




FIG. 10. The change of the translational energy of the par- 
ticle i due to sliding (AEl Up in the text). The solid line shows 
the result applicable to symmetric collisions (e n = 0.9 and 
C = 0.13). 

Figure 10 shows the results for AE l slip . We chose to 
show the total energy change, instead of the relative one, 
in order to be able to address the case of initially sta- 
tionary particle, characterized by Eq = 0. The solid line 
shows the result for the symmetric case, v p = 0. From 
Fig. 10 we observe that the loss of energy of the particle 
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strongly depends on v p , i.e. on the degree of the asym- 
metry of the collision. In particular, we observe that 
AE l slip could attain negative values, meaning that the 
particle increases its translational kinetic energy due to 
sliding. In order to illustrate this rather counter-intuitive 
point, let us consider for a moment completely asymmet- 
ric collision, characterized by v° = v°i, v° = 0. Using 
Eqs. ( |60| , [A?]), the change of the energy of the particle 
i (initially stationary particle), due to sliding, easily fol- 
lows 



AE 



slip 



\{l + e n )-C 
(1 



nB, 2 y 



. 1 mR 2 
(l + e„) ( - + — — 



-C 



(v [> 



(62) 



Since C = O(e) <C 1, AE l slip is positive, meaning that 
the particle i loses its translational energy due to sliding 
after the collision. On the other hand, the change of the 
energy of the particle j (impact particle), due to sliding, 
is given by 
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The negative sign implies that the particle j gains trans- 
lational energy by sliding. The interpretation of this re- 
sult is simple, in particular in completely elastic limit, 
e„ — > 1 (also C — > 0). Since the collision is elastic, the 
translational velocity of the impact particle j vanishes 
immediately after the collision with the stationary par- 
ticle i. But, the particle j still has the angular velocity, 
ftj, which is (in the elastic limit) equal to its initial angu- 
lar velocity. Consequently, the particle j has the sliding 
velocity, which is, immediately after the collision, equal 
to the negative of its initial translational velocity. The 
sliding acceleration resulting from this sliding velocity in- 
duces the motion of the particle in its initial, i, direction. 
The result is that the translational energy of the parti- 
cle j is being increased by the action of the friction force 
between the particle and the substrate after the collision. 

Still considering completely asymmetric case, it is of 
interest to compute the net energy loss of the system of 
two particles, AE l J ip = AE\ lip + AE 3 slip . By combining 
Eqs. (pi 1631) , we obtain 



asyrnm |(1 + e„) — C 



mR 2 ' 



(1 



mR 



(1 - e«) 1 + 



mR 2 



(v°) 2 



(64) 



relative loss of energy by dividing with the total initial 
translational kinetic energy (scaled with reduced mass), 
AE l s f ip = A.E z s f ip /(v ) 2 . I n the completely elastic case, 
the result for the relative loss of energy is given by 
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Following the same approach, the relative loss of energy 
of the system of two particles undergoing a symmetric 
elastic collision (scaled with the total initial energy) is 
given by (using Eqs. ([|| |f)) 
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Comparing Eqs. ( p5| , 66), we see that the particles lose 
twice as much energy due to sliding in symmetric, com- 
pared to completely asymmetric elastic collision. The 
intuitive understanding of this result follows by realiz- 
ing that the sliding velocities of the particles at the end 
of a symmetric collision, scaled by the initial velocities, 
are larger in the symmetric, compared to the completely 
asymmetric case (viz. Eq. (j57|)). The consequence is that 
the particles that have undergone a symmetric collision 
slide longer and lose more translational energy. When 
C^O, the loss of energy due to sliding in an inelastic col- 
lision is even smaller, since the particle-particle interac- 
tion during the collision decreases the angular velocities 
and, consequently, the sliding velocities of the particles 
after the collision. 

Figure 11 shows the change of momentum due to slid- 



ing, defined as Ap 



slip 



(v 



/ 



i, so that it measures 



the change of the translational velocity of particle i (in 
the i direction), after the collision. Clearly, Ap l slip de- 
pends very weakly on the degree of the asymmetry. For 
completely elastic collisions, Ap l slip depends only on the 
relative initial velocity of the particles, and it is given by 
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(67) 




The net change of the translational energy is positive, 
as expected, so that the system is losing translational 
kinetic energy. As in the symmetric case, we obtain the 



FIG. 11. The change of the linear momentum of the parti- 
cle i due to sliding. The solid line shows the result applicable 
to symmetric collisions (e„ = 0.9 and C = 0.13). 
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Effective coefficient of restitution. Let us define t m 
as the time, measured from the end of a collision, at 
which neither of the particles slides anymore, so that 
t m = max(tf, tj), where t s k (k = i,j) is the sliding time 
of a particle, given by Eq. (p9). We define the effective 
coefficient of restitution, ef/\as the ratio of the transla- 
tional velocities of the particles at the time t m , and their 
initial velocities. Using the translational velocity of the 
particle i, given by Eq. (|6C]), and the analogous equation 
for the particle j, we obtain 

|Vi-Vj|° 1 + 22^ 

Remarkably enough, this result involves only the "real" 
coefficient of restitution and the geometric properties of 
the particles. For solid spheres, the difference between 
the usual coefficient of restitution and the effective one is 
huge; for e„ = 0.9, we obtain e%ff = 0.36. This value is 
smaller than the range reported in [|TJ , but very close to 
our experimental results for steel particles on aluminum 
substrate |22}| . Slight imperfections from the spherical 
shape in experiments, non-central collisions, and/or the 
fact that the static friction between the particles has been 
neglected in our calculations, might be the reason for this 
discrepancy. 

General remarks. While more precise analysis and ma- 
terial parameters could be used in order to more precisely 
model experiments, we consider that the main results and 
observations given in this section are model- independent. 
In particular, the observation that the sliding is likely to 
occur as a consequence of most of the collisions does not 
depend on the details of the model. Of course, the re- 
sults would be modified in the case of more complicated 
(two dimensional) geometry of the collisions. Still, the 
particular geometry of a collision enters into our results 
for the energy and momentum change only through the 
observation that the frictional interaction of the colliding 
particles with the substrate can be ignored during a colli- 
sion. Since, for the system that we consider in this work, 
the collision forces are generally much stronger than the 
friction forces resulting from particle-substrate interac- 
tion, we do not expect this observation to be modified 
for more complicated collisions. We do note that more 
realistic model for the particle-particle interactions (e.g. 
by including static friction) would introduce modifica- 
tions in the expression for the final angular velocity of 
the particles, Eq. (|56|). 

In the experiments |^l|,^2| it is observed that some of 
the particles travel for long distances without colliding. 
Especially in this situation, it is important to include the 
effect of rolling friction, which we have ignored in this 
section. As long as a particle slides, the effect of rolling 
friction could be safely neglected, since the coefficient 
of rolling friction is much smaller than the coefficient of 
kinematic sliding friction. 



V. CONCLUSION 

The most important observation made in this work is 
that sliding leads to a considerable modification of the 
translational kinetic energy and linear momentum of the 
particles, even in the limit of completely elastic collisions. 
Based on this observation we give the result for the "ef- 
fective" coefficient of restitution, valid for dilute systems, 
where the mean free path of the particles in between of 
the collisions is much longer than the sliding distance. 
For more dense systems, we conjecture that this "effec- 
tive" coefficient of restitution strongly depends on the 
local density and velocity of the particles. 

The model that we present is to be used in molecular 
dynamics (MD) type simulations of externally driven 
system of a set of particles interacting on a horizontally 
oscillated surface. In particular, we have prepared the 
grounds for detailed modeling of the system of two kinds 
of particles, which are characterized by different rolling 
properties. In (22) is shown that strong segregation can 
be achieved. Preliminary MD results, based on the model 
formulated in this paper, show that the realistic modeling 
of the particle-particle and particle-substrate interactions 
are needed in order to fully understand this effect. 

Further, since experiment is the ultimate test for every 
theory, it would be of considerable importance to extend 
the previous work [T^-pO]] on formulating continuum the- 
ory for "2D granular gas" . Using the model presented 
here should allow for precise comparison between exper- 
imental and theoretical results. Possible formulation of 
realistic continuum, hydrodynamic theory applicable to 
this seemingly simple system would be an important step 
towards better understanding of granular materials. 
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APPENDIX A: LINEAR MODEL FOR THE 
NORMAL FORCE BETWEEN PARTICLES 

Let us analyze a simple situation, a central collision of 
two identical particles, i and j, moving with the veloci- 
ties, and v9, in the i direction only. Here we ignore 
the interaction of the particles with the substrate; the im- 
portance of this interaction is discussed in Appendix D. 
Using this assumption, the normal force, given by Eq. (0), 
is the only force acting on the particle i in the normal 
direction. By combining the equations of motion for the 
particles i and j, we obtain that the compression depth 
x = (d — Tij)/d satisfies the following equation 

x + 7at± + uj^x — , (Al) 
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where is the damping coefficient in the normal direc- 
tion, and ojq = y/2k/m. We limit our discussion to the 
case of low damping, so that e = jn/^o "C 1. 

This equation is subject to the following initial condi- 
tions: x(t = 0) = 0, x(t = 0) = v® el /d. The relative ve- 
locity of the particles at t = is given by v® el = v° — v° | ; 
for a symmetric collision, v® el — 2v°. The solution is 




V w o a - (¥) 
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The duration of the collision, t co i, now follows from the 
requirement x(t = t co i) = 0, thus 
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In what follows, we also need the time of maximum com- 
pression, t max . From the condition x(t = t max ) = 0, we 
obtain 
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Expanding to 0(e), it follows 
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(A5) 



So, the damping manifests itself in a slight asymmetry 
of the collision, since t max < t co i/2. The maximum com- 
pression, x max = x(t = tmax), follows using Eq. flA^). 
It is given by Eq. ([59|) for inelastic collisions, and by 
Eq. ( |32| ) for elastic ones. 

We define the coefficient of restitution as the ratio of 
the final velocities of the particles relative to their initial 
velocities, i.e. e„ = |vj — Vj|^/|vj — Vj\°. It follows 



APPENDIX B: NONLINEAR MODELS FOR THE 
NORMAL FORCE BETWEEN PARTICLES 

The linear model, presented in the previous section, is 
the simplest approximation for the collision interaction 
between particles. Nonlinear terms, resulting from the 
final area of contact and other effects 
should be included in order to model the interaction be- 
tween particles more realistically. Still, since we are con- 
cerned with the conditions where the maximum compres- 
sion depth is small, the linear approximation is a reason- 
able one. We use the nonlinear model, outlined below, in 
order to connect the values of the parameters, in partic- 
ular the collision time, t co i, with the material properties 
of the particles. 

The general, commonly used equation is [p0[ 



x H x 1 x + 
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where 77 and E are the material constants. The choice 
7 = 0, (3 = 0.5 leads to the Hertz model. The analysis of 
this equation gives an expression for t co h that can then be 
used to determine the appropriate force constant in the 
linear model, fc, and the damping coefficient, 7^. The 
result for the collision time is l|30| 
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For the Hertz model, 7(0.5) = 2.94. The parameter E is 
given by Y/(3(l — ct 2 )), where Y is the Young modulus, 
and a is the Poisson ratio. We use Y = 2.06 x 10 12 
dyn/cm 2 , and a = 0.28. For steel spheres with diameter 
d = 4 mm, and impact velocity, vq = 10 cm/s, f co ; « 
2.55 x 10~ 5 sec; for v = 100 cm/s, t col « 1.61 x 10" 5 

— 1/5 

sec. We note that the model predicts t co i ~ w , and 
tcoi ~ d. The parameters that enter the linear model 
can now be calculated, using u>q = 7r/t co ;(l + (9(e 2 )), and 
7at = ~2/t co i ln(e„). 
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In the limit of low damping, e n is close to 1; typically 
we use e„ = 0.9, appropriate for steel particles [^l] p3[ . 
Using Eq. (|A^), we obtain e w -2/7rln(e„) « 0.07. 

The final velocity of the particle i (at the end of the 
collision) follows from the requirement that the total lin- 
ear momentum is conserved in the center of mass frame. 
It is given by 



vf = ^[v° + v°- e „(v° "° 



(A7) 



For a symmetric collision, this results simplifies to 



-e v° 



APPENDIX C: SLIDING DURING COLLISIONS 
1. Sliding during a symmetric collision 

In Appendices [A| and |b] we obtained the results govern- 
ing dynamics of particle collisions, ignoring the interac- 
tion with the substrate. Here we show that the colliding 
particles slide through most of a typical collision. The 
additional material constants that are involved are the 
coefficients of static and kinematic friction between the 
considered particles and the substrate, /i s and fi^. In our 
estimates, we use fj, s = 0.5 and /ij, = 0.1. 

The condition for sliding, Eq. (|25|), applied to the sim- 
ple situation outlined in section [IV A| , gives that sliding 
occurs when |F^| > (1 + mR 2 / 1) n s mg. In terms of the 
compression depth and velocity, this condition is 
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du>QX + drf 
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X > 1 + 
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We note that the left hand side of this equation is always 
non-negative, since ¥ N is always repulsive (at the very 
end of a collision, when x <C 1, x < 0, F C N is set to 0). 
In the limit — > 0, we obtain that sliding occurs when 
x > x° mm , where x° min = (1 + mR 2 / I)[i s g/Jduj%). Using 
the result for the compression depth, Eq. (A2), we obtain 
the time at which sliding starts, ^ in , measured from the 
beginning of the collision, 



sin(oj t° mm ) = 1 + 
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For the initial velocities, v® el , satisfying v® el ^> v b , where 
v b = (1+mR 2 /I)/j s g/uj , it follows that sin(u/o^ in ) <C 1. 
For our set of parameters, and assuming solid spheres, 
v b ss 10~ 2 cm/s. Therefore, this condition is satisfied for 
most of the collisions. Assuming this, we obtain 
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Exploiting the symmetry of an elastic collision, we con- 
clude that the sliding condition is satisfied for f*U n < t < 

tool ~ train' 

Next we go to the limit of small, but finite damping, 
and assume that the condition uotmm <SC 1 is still valid, 
where t m i n is now the time when sliding occurs for ^ 
0. Using jNt m in < C LUot m i n , we Taylor-expand x and x 
(given by Eq. (A2)) at t = t m i n , and keep only the first 
order terms in small quantities toot m i n , jNt m in- I n this 
limit. 
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The sliding condition, Eq. (CI), gives the time when slid- 
ing occurs, for an inelastic collision 



t n 



1 



mR z 



UJ 



(C6) 



that the friction term could be relevant only in the limit 
e n — > 1, since v c diverges in this limit. Consequently, it 
follows that t m in — » 0, so that the sliding starts immedi- 
ately at the beginning of an inelastic symmetric colli sion . 
Since t m in — > 0, the expansion used to obtain Eq. (C5) 
is consistent. 



2. Sliding during an asymmetric collision 

By combining Eqs. ( [i^ , j54|) , we obtain the condition 
for sliding during an asymmetric collision 
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Using Eqs. (||, [35]) for F C N i and F C R -, respectively, we ob- 
tain (in terms of the compression depth, see Appendix ^) 
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where v* el is given by Eq. (Q). From the first part of 
this Appendix, we already know that the first term on 
the right hand side is negligible. The term inside the 
square brackets is positive for solid sphe res, and /i s = 0.5. 
For large x's, the condition, Eq. flC§|), is always satis- 
fied, since duj 2 x is the dominant term. So, we need to 
explore only the beginning and end of a collision. If 
sign{—v* el ) < 0, the sliding condition is always satis- 
fied; so that the slower particle always slides. When 
sign{—v* el ) > 0, we concentrate on the very beginning 
of the collision, and obtain the condition 
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Since typically 7s = 7at/2, this condition is satisfied, 
assuming |v° + v°| « |v? — v°|. 

We conclude that the particles entering an asymmetric 
collision slide during the whole course of the collision, 
except possibly in the case |v° — v9| <C |v? 
not consider this case here. 



Vj u |. We do 



We note that there are two factors that contribute to 
tmin'- the frictional interaction with the subst rate gives 
the first term on the right hand side of Eq. (C6), and 
the damping that occurs during a collision gives the sec- 
ond one. For the initial velocities, satisfying v° el ^> v c = 
g^s/lN, the contribution from the damping is the im- 
portant one. Using the expression for 77V given in Ap- 
pendix [b| we obtain v c as 0.05 cm/s (for e„ = 0.9). This 
velocity is smaller than the usual initial velocities con- 
sidered in this work. Assuming v® el » v c , we conclude 



APPENDIX D: MODIFICATION OF COLLISION 
DYNAMICS DUE TO THE INTERACTION 
WITH THE SUBSTRATE 

Here we estimate the importance of the interaction be- 
tween the colliding particles and the substrate during a 
collision. In particular, we estimate under what condi- 
tions the interaction with the substrate significantly mod- 
ifies the results for the compression depth and the dura- 
tion of a collision. We use the linear model outlined in 
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Appendix [X[ and concentrate on the case of the particles 
moving on a horizontal static substrate. 

In Appendix ^| it is shown that, assuming typical ex- 
perimental conditions, the colliding particles slide rela- 
tive to the substrate during most of a collision. For sim- 
plicity, here we concentrate on a symmetric collision, and 
further assume that the condition for sliding is satisfied 
throughout the collision, so that the friction force attains 
its maximum allowed value, given by Eq. ( |l8| ) . By using 
this approximation, we slightly overestimate the influ- 
ence of the friction with the substrate on the dynamics 
of a collision. 

From Fig. 4 we observe that the friction force, f , acts 
in the direction opposite to the normal collision force, 
F%. Including f in the Newton equations of motion for 
the particles i and j, we obtain the modified equation for 
the compression depth 
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which simplifies to Eq. (JA1J) if the particle-substrate in- 
teraction is ignored. 

Using the initial conditions as in Appendix |a|, we ob- 
tain the solution 
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where Xf = /ifcg/(cL->o). 

Collision time. For simplicity, we concentrate on 
the case of zero damping (■jn = 0), and calculate the 
change of the duration of the collision due to the particle- 
substrate interaction. Let us assume that the change of 
the collision time is small, and write t' col — i co z + T, where 
t C oi = ?r/a;o is the collision time if there is no interaction 
with the substrate, and r t co i- Using the condition 
x(t = t' col ) = 0, and expanding the compression depth, 



given by Eq. (D2), to the first order in the small quan- 
tity tloq, we obtain that r = 2x fd/v® el . So, the relative 
change of the collision time due to the interaction with 
the substrate is given by 
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For v® el 3> v a , where v a w ^i k g/ui , the change of the 
collision time is small. Using the parameters given in 
Appendices |b] and |c], we estimate v a rj 10 -3 cm/s. So, 
for most of the experimentally realizable conditions, the 
duration of a collision is just very weakly influenced by 
the particle-substrate interaction. We assume v® pl ^> v a , 
so that tuiq <C 1, and the expansion of Eq. (D2) is con- 
sistent. 

Maximum compression depth. Following the same ap- 
proach, we estimate the modification of the maximum 



compression achieved during a collision, due to the inter- 
action with the substrate. Working in the limit of zero 
damping, and assuming a small modification of the time, 
tmax, when the maximum compression, x' max , is reached, 
we obtain x' max « x f + v rei/ (du>o). Comparing this result 
with the result for the compression depth calculated pre- 
viously, given by the elastic limit of Eq. (p9), we obtain 



dujt 



x f 



Vk9 



(D4) 



Similarly to the analysis of the collision time, we observe 
that for v® el 3> v a , the maximum compression depth is 
very weakly influenced by the particle-substrate interac- 
tion. 

We conclude that for most of collisions occurring in ex- 
periments, the interaction with the substrate just slightly 
modifies the compression depth and the duration of a 
collision. These small modifications are ignored in the 
subsequent analysis. 



APPENDIX E: ROTATIONS OF THE PARTICLES 
DURING A COLLISION 

1. Rotations during symmetric collisions 

During symmetric collisions, the rotational motion of 
the particles is influenced only by the friction force be- 
tween the particles and the substrate. Here we consider 
only elastic collisions, since in Appendix ^ it is shown 
that the particles entering an inelastic collisions start 
sliding immediately, so that the angular acceleration is 
constant during the whole cours e of co llision, simplify- 
ing the calculations (see section I V A 2 ) . Since there is 
no possibility of confusion, we use scalar notation, with 
the sign convention that +sign corresponds to the forces 
acting in +i direction, and to the angular motion in +j 
direction (the coordinate axes are as shown in Fig. 4). 

At the very beginning of an cla stic collision, for < 
t < t a min , (t° min is given by Eq. (|C3|)), the colliding parti- 
cles do not slide. During this time interval, the angular 
acceleration of the particle i, which initially moves in — i 
direction, is given by 



R f R N, i 

T Jl 



I 1 + inf 



R 1 



7nR 

-l—.v^lt , (El) 



where x w tv® el /d, and Eqs. (||, |3l]) have been used. In- 
tegration yields 

W = tLn)=V + l^(l + ^)^, (E2) 
2 I V I J Vrel^O 

and fi? = -v°/R. For t° min < t < t co i - t° nml the sliding 
condition, |fj| = /i s |F^ J, is satisfied, so that the angular 
acceleration reaches its maximum (constant) value 
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a, 



mR 



(E3) 



For t co i > t > tool — tmim the sliding condition is not 
satisfied anymore, but the particle is sliding already, so 
that (li is still given by Eq. (E3). The angular velocity 
of the particle i at the end of collision is 



as it follows from Eq. (|51|), where is being neglected. 
In performing the integration, it appears that we have 
to consider separately two regions: < t < t crlt , dur- 
ing which FJj varies, and t > t crlt , during which F^ is 
constant. The final angular velocity of the particle i is 
formally given by 



Cl 



. mR 



t. 



(E4) 



Combining Eqs. (|e|, HJ) , we obtain the final result, given 
by Eq. A 



2. Rotations during asymmetric collisions 



a. About tangential force 



Here we estimate under what conditions, F C R , given by 
Eq. (0), reaches its m axim um allowed value, i/fc|Fjy-|. As 
mentioned in section [VB, here we ignore the frictional 
interaction of the particles with the substrate during a 
collision. For simplicity, we also neglect the damping 
in the normal directions, so that |F^| = mdu^x (see 



Appendix [A]) . Next, we note that the relative velocity of 
the point of contact satisfies v z el (t = 0) > v* el Ct > 0), 
since F^, always decreases v* el (given by Eq. ([?])). In 
what follows, we use v* el (t > 0) = v z rel {t = 0), and give 
the upper limit of the first term entering the definition 
ofFfc. 

Let us first concentrate on large compression depths, 
v°|/(dw ) (see Appendix 0). This 



compression is reached at t = t r< 



7r/(2w ). We use 



l \P ^max) 



rel 



(t = 0) = |v° + v^|, and obtain that 



F C R reaches its maximum allowed value if (see Eq. 



7s 
2 ' 



(E5) 



Since 75 /loq <C 1, this condition is never satisfied for 
v k = 0{l), and |v?+v?|«|v?-v?|. 



For small x's, let us assume again v* el (t > 0) 



rel 



(* = 



0). From Eq. (|6|) it follows that F^ reaches its cutoff 
value when x < x crlt , where 



n/ = n? 



Cljdt 



Cljdt 



(E8) 



This result can be simplified by realizing that |F^| = 
O(e). It follows that \Cli\ = O(e), so that the contribu- 
tion of the second term on the right hand side of Eq. (ES ) 
is proportional to \Cli\t crlt — 0(e 2 ). For consistency rea- 
sons we neglect this correction, and ignore the fact that 
|Fjj| could reach Coulomb cutoff at the very beginning 
and end of a collision. This estimate is not valid for 



-vv| < ivy 



-Vj\, so when the particles initially move 



with almost the same velocities. As already mentioned 
in Appendix [c], we do not consider this case here. 



b. The angular velocity of the particles during an 
asymmetric collision 

Using Eqs. ( [5l| , j55[ ), and neglecting the particle- 
substrate interaction during a collision, we obtain the 
angular acceleration of the particle i, 



Cl; 



mR 2 7s 
7 2 ' 



:a + «,■)■]] l 



(E9) 



and Clj = Cli. Recalling that Cli and Clj are always in 
the opposite direction from Cli + Clj , we obtain a simple 
system of coupled ordinary differential equations 



Cl, = -C'iCli + Clj) 



Jj > 



CI, 



(E10) 
(Ell) 



where C = mi? 2 7 S /(2J). We define Cl + = Cli + Clj, so 
that Cl + = — 2CCI+, with the solution 



Cl + (t) = Cl+(t = 0) exp(-2C"i) 



(E12) 



7s 



2v k uj Q |v J 



0(e) 



(E6) 



Using x w t\v® — Vj\/d (valid for x <C x max ), we obtain 
that the condition x < x crlt is satisfied for t < t crlt , 
where 



7s 



7TW0 |V? 



= 0(e) . 



(E7) 



In order to calculate the angular velocity of the parti- 
cle i at the end of a collision, Cl{, we have to integrate 
the angular acceleration, Cli, during the course of a col- 
lision. The angular acceleration is proportional to F C R , 



At t = tcoi, Cl + (t = tcoi) = Cl + (t = 0)exp(-2C), where 
C = C t co i = 0(e). Recalling that the changes of Cli and 
Clj are the same, so that Clk(t = t co i) = Clk(t = Q) + ACl, 
(k = i, j), the change of the angular velocities is given by 

ACl = \ [(fl° + f2°)(exp(-2C") - 1)] « -(n? + Cl°)C , 

(E13) 

correct to the first order in e. For 7s = 7jv/2, and the 
parameters as in Appendix [§], C = nmR 2 e j (47) w 0.13. 
The final angular velocity of the particle i is now given 
by Eq. ©. 
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APPENDIX F: JUMP CONDITION FOR 
ASYMMETRIC PARTICLE COLLISIONS 

Throughout this work, we have assumed that the par- 
ticles are bound to move on the surface of the substrate. 
Here we explore the validity of this assumption. The 
required condition for a particle to be bound to the sub- 
strate is that the normal force |Fjy|, given by Eq. (54), 
is nonzero. We immediately observe that only a particle 
colliding with a slower particle (so that sign(—v z el ) < 0, 
see Eq. (Q)) experiences a force in +k direction, due to a 
collision. Let us concentrate on this situation. Using the 
value of |Pjj| at t = 0, we obtain that a particle detaches 
from the substrate if 



fn 



IS i 



9) >0 



(Fl) 



where Eqs. (0, |54|, [55]) have been used. It follows that, 
during the collisions distinguished by \v z el \ > v d = 
2g/js, the faster particle detaches from the substrate. 
Using the values of the parameters as in Appendix 
and 7s = 7at/2, we obtain v d w 0.5 cm/s. Correspond- 
ingly, this effect takes place during most of the asymmet- 
ric collisions occurring in typical experiments (2^^^] . By 
relating the impulse of the force / ne t transferred to a par- 
ticle while the collision is taking place, with the change 
of the momentum of the particle in the k direction, we 
obtain the estimate for the initial velocity of the particle 
in the k direction 



7T 

9) — 

UJQ 



(F2) 



The maximum height above the substrate which the par- 
ticle reaches is h z = l/2(v z ) 2 /g, and the time spent 
without contact with the substrate is t z = 2v z / g. Let 
us assume a completely asymmetric collision, so that 
|vf | = v°, |v°| = 0, and v re i — v z el — v° . Using the 
parameters from Appendix [b], for v° = 10 cm/s, we ob- 
tain v z w 0.5 cm/s, h z « 1.3 x 10 -4 cm, and t z w 10~ 2 
sec. Since the maximum height is much smaller than 
the diameter of the particles, this detachment introduces 
negligible corrections to the dynamics of the particle col- 
lisions in x — y plane. Further, even though t z 3> t co i, so 
that the particle is not in the contact with the substrate 
during the time which is much longer than the duration 
of the collision, t z is still much smaller then the sliding 
time scale, specified by Eq. (|59|). So, our results for slid- 
ing of the particles after a collision are not significantly 
modified due to the detachment effect. 
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